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Abstract 

O , 

The paper studies a particular class of analytic solutions for the Gen- 
eralized Ohm's Law, approached by means of the so called formal powers 
of the Pseudoanalytic Function Theory. The reader will find a description 

S of the electrical current distributions inside bounded domains, within in- 
homogeneous media, and their corresponding electric potentials near the 
boundary. Finally, it is described a technique for approaching separable- 
variables conductivity functions, a requisite when applying the construc- 
^ tive methods posed in this work. 

o 

00 1 Introduction 

<N 

y— | The study of the generalized Ohm's Law 

div (crgradu) = 0, (1) 

where a denotes the electrical conductivity function and u is the electric poten- 
• i-h tial, is the base for well understanding a wide class of problems in Electromag- 

netic Theory, just as the Electrical Impedance Tomography, the name given in 
medical imaging for the inverse problem posed by Calderon [3] in 1980. Yet, 
for many decades, the mathematical complexity of (JlJ imposed so difficult chal- 
lenges to the researchers, that the structure of its general solution in analytic 
form remained unknown. It was until 2006 that K. Astala and L. Palvarinta 
IT] discovered that the two-dimensional case of was closely related with a 
Vekua equation pj . , and in 2007 V. Kravchenko et al. i8. , based upon elements 
of Pseudoanalytic Function Theory [5] , achieved to pose what could be consid- 
ered the first general solution in analytic form of for the plain, when the 
conductivity function a belongs to a special class of functions. 

Virtually, these two discoveries were the departure point for developing a 
completely new theory for the generalized Ohm's Law, mainly because they 
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allowed to research a wide sort of electromagnetic phenomena that had remained 
out of range for the mathematical tools known previously. 

In this paper we discuss the possibility of considering the formal powers [2] , 
as a new scope for analyzing the electrical current distributions inside inhomo- 
geneous media in bounded domains, since it is bias their linear combination that 
we can approach the general solution for the two-dimensional case of pg, 
and they are also useful for constructing an infinite set of analytic solutions for 
its three-dimensional case [12) . 

Starting with the elements of Pseudoanalytic Function Theory, and of Quater- 
nionic Analysis, we briefly expose an idea for rewriting the three-dimensional 
case of ([I]) in a quaternionic equation, in order to pose the structure of its general 
solution by means of a generalization of the Bers generating pair, in Complex 
Analysis. 

Eventually, we focus our attention on the plane by considering one example 
in which a spacial variable is fixed, and the conductivity a adopts an exponential 
form. Then we trace the electrical current density patches, obtained from the 
solutions of the Vekua equation in terms of Taylor series in formal powers, and 
we show that, from an adequate point of view, the observed patches might keep 
a sort of regular dynamics when flowing through this inhomogeneous medium, 
once they are compared to those traced for an homogeneous case. 

The work closes with a basic idea for approaching separable-variables con- 
ductivity functions, since this is the central requirement if we are to apply the 
exposed mathematical ideas for approaching solutions of ([I]), and to analyze 
their meaning in terms of electrical current flows. 

2 Preliminaries 

2.1 Elements of pseudoanalytic functions 

Following L. Bers [2], a pair of complex- valued functions F = Re (F) + ilm (F) 
and G = Re (G) + ilm (G) will be called a generating pair if the following 
condition holds: 



Here i denotes the standard imaginary unit i 2 = — 1, whereas F represents the 
complex conjugation of F : F — Re (F) — ilm (F) . Thus, any complex-valued 
function W can be represented as the linear combination of the generating pair 



where <f> and ip are both real-valued functions. Hence the derivative in the sense 
of Bers, or [F, G)- derivative, of a complex- valued function W will be defined as 




(2) 



(F,G): 



W = <pF + 



(3) 



d (FjG) W = (d z (t>)F+(d^)G, 



(4) 



where d z 




and it will exist if and only if 



(0^) F + (fcV) G = 0, 



(5) 
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where c\ — + i ■ Notice even the operators d z and c\ are usually defined 
including the coefficient | , it will result somehow more convenient for this paper 
to work without it. 

By introducing the notations 

. _ Fd z G~G8 z F _ Fc\G-Gb\F 

A{F > G) - FG-FG ' a ^ G) ~ FG-FG ' (6) 
Fd z G - Gd z F _ Fc\G - Gc\F 

B {F>G) - F g_ Tg . b (F, G) - F g_f G i 

the (F, G)-derivative of Q can be written as 

9 ( f,g)W = A (F>G) W + B {F , G) W, (7) 
and the condition ^ will turn into 

d w W-a (FtG) W-b lP , a) W = 0. (8) 

The last differential equation is known as the Vekua equation [15j , and it will 
play a very important roll in our further discussions. The functions defined in 
([6| are known as the characteristic coefficients of the generating pair (F, G), and 
every complex- valued function W satisfying Q will be referred as an (F, G) - 
pseudoanalytic function. 

The following statements were originally posed in [2]. Another authors will 
be cited explicitly. 

Remark 1 The functions conforming the generating pair (F,G) are (F,G)- 
pseudoanalytic, and their (F,G) -derivatives are d( F _a)F — dr F> Q\G = 0. 

Remark 2 Let p be a non-vanishing function inside some domain CI. The pair 
of functions F = p and G = ~ satisfy the condition so they constitute a 
Bers generating pair, and their characteristic coefficients are 

A (F,G) = a (F,G) = 0, (9) 
r> 9zP , C\p 

B{ F ,G) - — , b {F ,G) - — • 

Theorem 1 Jfjj] Letp be a non-vanishing function within some domain CI. The 
pair of real-valued functions <f> and ip will be solutions of the system 

— 6-— — ib —<{)-- — —ih (10) 
dx p 2 dy ' dy p 2 dx 

if and only if W = p<p + ^iJj is solution of the Vekua equation 

c\W - —W = 0. 



A pair of functions <fi and ip satisfying (10) is called a p-analytic system 
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Definition 1 Let (Fq,Gq) and (F±,Gi) be two generating pairs, and let their 
characteristic coefficients satisfy 

a (F ,G ) = a (fi,Gi) and B (F ,G ) = -&(Gi,Ji)- (H) 

T/ie generating pair (i<\, Gi) wzH be called a successor pair of (F , Gq) , as well 
as the pair (F ,Gq) will be called a predecessor of (Fi,G\). 

Theorem 2 Let the complex-valued function W be (Fq, Go) -pseudoanalytic, and 
let the generating pair (F\,G\) be a successor pair of (F ,Gq). Hence, the 
(Fq, Go) -derivative ofW will be an (i*\,Gi) -pseudoanalytic function. 

Definition 2 Let the pairs of functions belonging to the set 

{(F n ,G n ):n = 0,±l,±2,...} (12) 

be all generating pairs, and let every (F n+ i, G n +i) be a successor pair of {F n , G n ). 



Then the set (12) is called a generating sequence. If the generating pair (F, G) 



(Fa, Go), we say that (F,G) is embedded in the generating sequence (12) 



Definition 3 The generating pairs (F ni G n ) and (F^, G' n ) are called equivalent 
if they posses the same characteristic coefficients |#|). 

Definition 4 A generating sequence {(F n ,G n )} is called periodic, with period 
a > 0. if the generating pairs (F n ,G n ) and (F n+a ,G n + a ) are equivalent. 

Definition 5 Let the complex-valued function W be (F,G) -pseudoanalytic, and 
let {(F ni G n )} be a generating sequence in which the generating pair (F,G) is 
embedded. The higher derivatives in the sense of Bers of W will be expressed as 

W [0] = W; 
W^=d iFniGn) W^; n = 0,1,2,... 

L. Bers also introduced the notion of the (F, G)-integral for a complex- valued 
function W . The following statements describe its structure, the necessary 
conditions for its existence and some of its properties. 

Definition 6 Let (Fq,Gq) be a generating pair. Its adjoint pair (F *,Gq) will 
be defined according to the formulas 



p* o q* _ 2( JQ 

FqGq — FqGq FqGq — FqGq 



In particular, if p is a non-vanishing function inside a domain Q, the adjoint 
pair of Fq — p, Gq = ^ will have the form 

F o* = -ip, Gq^P^ 1 - 
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Definition 7 The (F , Go)-integral of a complex-valued function W is defined 
as 



where A is any rectifiable curve inside some domain ft, going from z until z, 
and it will exist iff W satisfies 



Theorem 3 The (F , Go) -derivative of an (F , Go)-pseudoanalytic function W 
will be (Fo,Go)-integrable. 

Remark 3 If W = <j)F n + ipG n is an (F n , G n )-pseudoanalytic function inside 
ft, its (F n ,G n ) -derivative will be (F n ,G n )-integrable. Indeed 



where zq is a fixed point. Moreover, by Remark 1 the (F n ,G n )- derivatives of 
F n and G n vanish identically, thus the last integral expression represents the 
(F n , G„)-antiderivative of d(F n ,G n )W (z) . 

The following paragraphs will expose some definitions and properties of the 
so called formal powers, whose physical implications will provide most material 
for this work. 

Definition 8 The formal power Zn^ {ao, zq; z) belonging to the generating pair 
(F n ,G n ), with exponent 0, complex coefficient ao, center at the fixed point zq 
and depending upon z = x + iy, is defined by the linear combination of F n and 
G n according to the expression 



where A and [i are real constants such that 

\F n (zo) + [iG n (z ) = ao- 
The formal powers with higher exponents are defined by the recursive formulas 






d (F n ,G n )W (z) d {FntGn) z = W(z)-(f> {z ) F n (z) - V (z ) G n (z) , 



Z n 0) (a ,z ;z 



) = XF n + t iG 




where m, n = 0, 1, 2, ... 



Theorem 4 The formal powers posses the following properties: 

1) Z n m ^ (a mi zo;z) is {F n ,G n )-pseudoanalytic. 

2) If a' and a" are real constants, then 



ZW (a + ia", z ; z) = a'Z™ (1, z a ; z) + a" ' (i, z ; z 
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for m, n = 0, 1, 2, ... 

3) When z — > zq we have that 

^« m) («m> ^o; a m (z - z ) m . 

Theorem 5 Let W be an [F 1 G)-pseudoanalytic function. Then, it accepts the 
expansion 

oo 

W= ^Z( m )(a m ,z ;z), (13) 

m=0 

where the absence of the subindex "n " in the formal powers (a m , zq; z) im- 
plies that all formal powers belong to the same generating pair. The coefficients 
a m will be given by the formulas 

_ (z ) 
am. ; • 



The expansion (13) is called Taylor series in formal powers ofW. 

Remark 4 Since every [F, G) -pseudoanalytic function W can be expressed by 
means of the Taylor series {13]) , this expansion is in fact an analytical repre- 
sentation for the general solution of the Vekua equation Q) } being the formal 
powers 

Z (m) (l,z ;z) and Z {m) (i, z ; z) , 
by virtue of Theorem 15, number (2), a base for the set of its solutions. 

2.1.1 An alternative path for introducing the concept of Bers gen- 
erating pair 

Following [9], let <j> be a real- valued function and let F be a complex- valued 
function. Consider the equality 

(ch-a-bC)(<f>F) = {ch<t>)F- (14) 

where a and b are complex- valued functions, and C denotes the complex conju- 
gation operator acting upon F as CF = F. The partial differential operator in 
the left side of the equation is clearly the one corresponding to a Vekua equa- 
tion. A simple calculation will show that ( 14 ) will be valid if and only if the 
complex-valued function F is a particular solution of 

c\F - aF - b~F = 0. (15) 

In the same way, let ip be a real-valued function and let G be a particular 
solution of (15 1. Thus the equality 

(Bg-a-bC) (V>G) = (c^)G 
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will hold. By adding this equation with (14 1, and introducing the notation 
W — cf>F + ipG, we will obtain 



dzW -aW-bW = (dz<p) F + (d^>) G. 

Hence 

dzW - aW - bW = (16) 

if and only if 

(dri) F + (dsffi) G = 0. 

If we require the functions F and G to fulfil the condition |2| , we will arrive to 
the very definition of an (F, G)-pseudoanalytic function W. Indeed, additional 
calculations will show that the functions a and b are precisely the characteristic 
coefficients Q-(f,g) an d &(F,G) defined in (|6|), thus (16) coincides with 



This alternative procedure will be useful for introducing the notion of Bers 
generating sets for the solutions of the three-dimensional quaternionic General- 
ized Ohm's Law. 



2.2 Elements of Quaternionic Analysis 

The algebra of real quaternions will be denoted by H(K) (see e.g. [7J). Every 
element belonging to this set will have the form q = Y^k=o1k e ki where q^ = 
qk (xi, X2, X3) ; k = 0,3 are in general real- valued functions depending upon 
three spacial variables, eo = 1, and e^; k = 1,2,3 are the standard quaternionic 
units, possessing the following properties of multiplication: 

e? = e2 = e§ = -l, (17) 
eie 2 e 3 = -1. 



It will be useful to introduce the auxiliary notation for an element q € H (M) 

Q = Qo + ~$, 

where clearly ~~q* = X)fc=o^ fcefe ' Thus q will be named the scalar part of the 
quaternion q, whereas q will be referred as the vectorial part of q. It is important 
to point out that the set of purely vectorial quaternionic functions such that 
{q = ~cf : q € H (R)} conforms an isomorphism with the set of three-dimensional 
Cartesian vectors R 3 . 



As it can be easily inferred from (17), the quaternionic product is not com- 



mutative, hence the multiplication by the right hand-side of the quaternion p by 
the quaternion q will be written as 



pq = M q p. 
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2.2.1 The Moisil-Theodoresco differential operator 

On the set of at least once-differentiable quaternionic-valued functions, it is 
defined the Moisil-Theodoresco partial differential operator, that was first in- 
troduced by Hamilton itself 

3 „ 

k=l K 

By means of the isomorphism remarked before, the operator D acts upon a 
function q £ H (R) according to the rule: 

Dq = gradgo — divq + rot "if, (18) 

where "grad", "div" and "rot" are the classical Cartesian operators, written 
using quaternionic notations (see e.g. [7]). 

2.2.2 Bers generating sets for the solutions of partial differential 
equations 

Let us consider the quaternionic differential equation 

DQ + Qp = 0, (19) 

where p,Q £ M (R) . Employing the idea exposed within the last paragraphs 
dedicated to the elements of pseudoanalytic functions, let Q be a particular 
solution of the equation ( 19 1 and let <p be a purely scalar function. A short 
calculation will show that the following equation holds 

D(<pQ) + {<pQ)p = (D<p)Q. 

Hence, if we own a set of four linearly independent Q^gi (R) ; k — 0, 3, every 
one solution of (19), we will be able to represent the general solution of (19 1 by 
means of the linear combination of Qk ■ 

Theorem 6 Let {Qk\k = Q C H(R) be a set of linearly independent solutions 
of the equation (lSty , thus its general solution can be written as 

3 



Q = ^2<p k Qk, 

fc=0 

where {<v9fc}fc =0 C K are scalar functions, all solutions of 

:>, 

J2( Dl Pk)Qk = 0. (20) 



fe=0 



We shall call {Qfc}fc =0 ^ e Bers generating set for the solutions of (19). 
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3 Study of the generalized Ohm's Law for a separable- 
variables conductivity function a 

As it was shown in [12] and [14], introducing the notations 

£ = -y/agmdu, d = - — ^— ; (21) 
the Ohm's Law ([I]) can be rewritten in quaternionic form: 

' D + M^) ~t = 0. (22) 



Hence, by virtue of Theorem 18, its general solution will have the form 

3 



£ = ^ fk^k, (23) 



k=l 



where |^fc| C H(R) is a set of linearly independent solutions of (22), and 
{Pk}k=x C K are scalar functions, all solutions of 

3 

£k = 0. (24) 

fc=i 

Specifically, when the conductivity c is a separable-variables function cr = 



Si (xi) S2 {x2) S3 (x 3 ) the Bers generating set (23) can be constructed explicitly 
|12j (it is worth of mention that a separable-variables conductivity function, 
in polar coordinates, was considered in the interesting work jjj, for studying 
the Generalized Ohm's Law by means of different mathematical methods). For 
example, let £ x = ex£x, where £\ 6 R. Substituting into |22[ ), we will obtain 
the partial differential system 

— £ ! + flifc = 0, —Ex - s' 2 £i = 0, —Ex - s' 3 £x = 0, 
where s' k = ^gf^Sfc; & = 1,2,3, for which a solution is 

£ _ p - f s' 1 dx 1 +f s' 2 dx 2 +f s' 3 dx 3 

Using the same idea, we will find out that the set of quaternionic-valued func- 
tions 



tx 


= ex£x 


— eie - I s' 1 dx 1 +J s' 2 dx 2 +J s' 3 dx 3 


(25) 


? a 


= e 2 £ 2 


= e 2 eJ s 'l dx i-f s' 2 dx 2 +J s' 3 dx 3 




*, 


= e 3 £ 3 


= e 3 eJ S 'i dxi +I s 2 dx 2-f s' 3 dx 3 ^ 
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conforms a Bers generating set for the solutions of (22) 



We must find the general solution of ( 20 1 if we are to pose the general solution 



of ( 22 ) , but it is not clear yet how to achieve this task. However, it is possible to 



construct an infinite set of solutions for ( 22 1 by means of the following procedure 



Suppose the scalar function tp^ vanishes identically inside the domain of 



interest f2. Then (20 1 will turn into 

(Dtpx) ~2i + (D<p 2 ) t 2 = 0, 



which will reach the system 

d 1 d 

"~2 ^ — 
jr 0x2 

d 



dxi 



dx- 



dx 



d 1 

o = ~2 Jt 

OX2 P OX\ 

-y>2 = 0; 



(26) 



where p = e~ J s i dx ^+f s 2 dx 2 _ The first pair of equations conforms precisely the 
p-analytic system introduced in Theorem 3, thus its corresponding Vekua 
equation will have the form 



P 



W = 0, 



(27) 



where W = p Vl + and t\ = ^ + 1^. 

Moreover, based upon a result posed by L. Bers itself in [2], and latter 
generalized by V. Kravchenko in 6J, slightly adapted for this work, we are able 
to find in analytic form the generating sequence that will allow us to build the 



formal powers for approaching the general solution of (27) in terms of Taylor 
series. 

Theorem 7 [21 [61 Let p be a non-vanishing separable-variables function inside 
a domain O, and let Fo = p, Gq = ^. Since p = p\ {x\)p2 (x 2 ) the generating 
pair (Fq,Gq) is embedded within a periodic generating sequence {(F m ,G m )}, 
with period 2, such that 



TP 

-*- n 



Pi (xi)P2 (^2) , G Tl 



when m is an even number, and 

Pi (gi) 
Pi {x 2 ) ' 

when m is odd. 



TP 



G r , 



Pi (Xl)p2 (X2) ' 
.P2 (X 2 ) 



Pi (xi) 



It is evident that identical procedures can be employed for the cases when 



the scalar function ipi vanishes identically in (20), and when does; leading to 



an infinite set of solutions for the three-dimensional quaternionic Generalized 



Ohm's Law ( 22 ) 



10 



3.1 Electrical current distributions in analytic form 

In order to pose an example for analyzing the behavior of the electrical cur- 
rent distributions, obtained bias the formal powers method; let us consider a 
conductivity function a with the form 

a = Sl ( Xl ) s 2 (x 2 ) s 3 (a*) - e 2 ^ +2 ^+ 2 ^, 
where a\,a 2 and 173 are all real constants, and let tp3 vanish identically inside 



Q,. In concordance with (26) the function p will have the form 



and the corresponding Vekua equation will be 

d^W- (02 -iax) W = 0, (28) 

where cV = -J^ — h i-^—. Following the Theorems and Definitions cited in the 
section of Preliminaries, let us consider the domain Q as the unitary circle, 
fixing the center of the formal powers at the origin of the plane Co = (it is 
remarkable that an equivalent case was considered in [1], where the authors 
performed numerical calculations in order to examine solutions for boundary 
value problems for elliptic differential operators). 

Since every formal power Z^™^ (a m , 0; £) > C = X2+IX1, is (F, G)-pseudoanalytic, 
it will be enough to focus our attention into the electrical currents emerging from 
the pair of functions Z^ m ~> (1,0; C) and Z^ (i,0;C) , = 0, 1,2, because any 
other electrical current will be necessarily a linear combination of them (see 
Theorem 15). 

Hence, let us consider the three first pairs of formal powers, corresponding 
to the generating pair Fq — e ~' TlXl+ ' J2X2 and Go = ie rTlXl ~ a ' 2X2 , in exact form: 

Z (0) (1,0; C) = e- aiXl+aaX2 , Z m (»,0;C) = ie* 1 * 1 "" 12 ; (29) 



Z (1) (1,0;C) = —e- TlXl smh(a 2 x 2 ) + —e-' T2X2 smh((j 1 x 1 ), (30) 

Z (1) (i,0;C) = - — e ff2!t2 sinh(CT 1 x 1 ) + — e a 1X1 sinh (cr 2 x 2 ) ; 
U\ CT 2 
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and 



Z <«(l,0;C)=(i + 



X2_ 

-1T2^ 2X2 sinh [cxtxx) - ^e^ 1 * 1 sinh (a 2 x 2 ) + 
(e aiXl sinh (<j 2 x 2 ) - e- a ^ sinh (a lXl )) + 

2aia 2 

+ - l — s i n h (aixi - a 2 x 2 ) , 

2<7 1 xi - 2a 2 x 2 \<T2 <?! J 

+ ^2 e"" 2 " 2 sinh iWi) + ^ aiXl sinh (ct 2 x 2 ) + 

1 2 

+ (e- CTlXl sinh (a 2 a; 2 ) - e^ 2 sinh (a lXl )) + 

Lo\<j 2 

+ - — [ — - — ] sinh(crij:i - d 2 a; 2 ) . 

2a 1 x 1 -2a 2 x 2 \a 2 o\) 



(31) 



Taking into account the notations introduced in (21 1, (23 1, (25) and (27) 
we have that the current density vectors j l - m ' will be given by 

y /aSip- 1 RjeZ<- ni 
j >! •' , = -a^rndu = \ ^a£ 2 phnZ^ ) . i*2) 

3.2 Electrical current patches within the unitary circle 

In order to pose a qualitative idea of the current patches inside the unitary 
circle, it will be convenient to observe first the traces provoked by the electrical 
currents flowing through an homogeneous medium, say a — I. Evidently, this 



implies a\ = a 2 = 0, and thus the Vekua equation (28) will turn into the well 
known Cauchy-Riemann equation 

d^Wh = 0, 

for which the standard Taylor series describe the general solution 
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-+(1) 

Figure 1: j h (1,0) 



By considering (32 1 for £°, £ and £ 2 , with coefficients 1 and i alternatively, we 
will obtain 



?i 0) (i,o)=l 



X-2 



?L 1} (i,o)= I X! 



?L 2) (i, o) 




?i 0) M) 



?L 2) M) 




(33) 



To trace the electrical current patches for (1,0) and j ^ (i, 0) is a triv- 
ial task, since they posses only one spatial component. Some of the patches 
corresponding to the rest of vectors are shown in the following figures. 

In these Figures, as for all corresponding to this section, the arrows point 
out the direction of the electrical current flows. Just such nearest to the center 
are omitted by virtue of the space. 

Let us focus our attention now onto the inhomogeneous case. Just as it 
happened before, the traces corresponding to the currents emerging from the 
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-1 

-1 -0.8 -0 6 -0.4 -0.2 0.2 04 0.6 0.8 1 





->(2) 

Figure 4: j h (i, 0) 



formal powers (1, 0; £) and Z( ) (i, 0; £) does not reach interesting diagrams. 

( e 2cr 2 x 2 +2<j 3 x 3 \ / 

o , t (0) M;zi,* 2 )= e 2 <™+ 2 "»*3 

(34) 

More illustrative examples arise when considering the current flows of (1, 0; £) 
and (i, 0; £) (for simplicity, next Figures are traced by fixing x% = 0): 

/ ±e a2X2+2u3X3 smh(a 2 x 2 ) \ 
~j (1) (h0;x ll x 2 )= ( Je^^+^sinh^X!) j , (35) 



h2 CT2 x 2 +2 CT3 X3 sinll ^ ia;i ) 

(i,0; Xl ,x 2 ) = ( ^e 2CTlXl+CT2X2+2<73a:3 sinh(a2a;2) 





The diagrams are drawn considering a\ — 3 and <t 2 = 1; an d by evaluating 
the current density vectors j into a fixed point, displacing the mark in the 
direction of the vector in a distance of 0.01% of the Cartesian norm in the 
subregions with higher conductivity, and 10% in such with lower. 

The Figures 7 and 8 show current patches when considering the remaining 
formal powers j ^ (1, 0; x\,x 2 ) and j ^ (i, 0; Xi,x 2 ). 

These short previews to the qualitative dynamics of the electrical current 
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Figure 5: j (1, 0; Xi, x 2 ) 




-1 -0.8 -0.6 -0.4 -0.2 0.2 0.4 0.6 0.8 1 



Figure 6: j (i, 0; Si, a; 2 ) 
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density vectors j ( m \ m = 1, 2; already allow us to point out some character- 
istics that seem to carry out useful patterns, when comparing them with the 
currents distributions for a constant conductivity function. 

For instance, Figure 5 keeps the dynamical behavior of the traces found 
in Figure 1, presenting the most significant variations into the positive X\- 
semiplane. Naturally, the magnitudes of the vector currents in both figures 
differ considerably when taking into account the conductivity values at every 
point into the domain S7. Still, a sort of dynamical relation might be suggested 
when adopting an appropriate point of view, as it will be posed in the next 
section, dedicated to the electric potentials at the boundary. 

It is also remarkable that already watching the traces of the very first formal 
powers, we can infer that the most important disparities between the current 
patches corresponding to the homogeneous case and such of the inhomogeneous, 
will take place on the diagrams belonging to the lower formal degrees. This is, 
to such formal powers whose superindex are among the smallest. 

This statement is based onto the location of what could be considered sink 
points and source points, following the terms employed in Complex Dynamical 
Systems (see e.g. [TU] )■ Of course, in a strict sense, not any sink or source point 
is found inside the unitary circle, neither will be outside of it, because all traces 
behave asymptotically. Nevertheless, the reader could easily infer the location 
of such regions by watching the places around the perimeter where the patches 
seem to get closer to each other. 

The bigger displacement of such regions takes place when appreciating the 
diagrams corresponding to the expressions with the smallest formal degrees. 
If the reader wishes to verify such behavior, perhaps standard computational 
languages for symbolic calculations can be useful in order to obtain in exact 
form some formal powers with higher degrees. Other case, it could be enough 
to trace again the posed diagrams considering, for example, o~% = o~% = 1. 

Thus, from this qualitative appreciation (at least for the inhomogeneous case 
just examined), we can observe a typical behavior of the formal powers when 
compared to the standard powers inside the unitary circle. We shall remember 
that this kind of behavior had been already remarked by L. Bers itself. Indeed, 
as it was mentioned before, professor Bers gave the complete proof for the case 
when z — > zo [2], opening the opportunity to study the behavior of the formal 
powers far away from their center. 

3.3 The electric potential at the boundary 

On the light of the last qualitative overview, a quantitative examination is 
in order. Let us pay our attention into one of the most important particular 
topics of the Generalized Ohm's Law, when studied into bounded domains: The 
behavior of the electric potential u at the boundary. 

It is easy to see that for approaching the electric potentials u^ m \xi, X2) 
corresponding to the current density vectors posed before, numerical methods 
are already needed for the cases corresponding to the formal degree 2. Let us 
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Figure 9: (1, 0; x u x 2 ) 
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Figure 10: (i, 0; x\, x 2 ) 



consider by now some details of the electric potentials that can be obtained in 
exact form. 

Ignoring the constant parts, the ele ctric potential corresponding to the cur- 



rent (1, 0; x\, X2) described in (34| is 



*(°)(l,0;xi,x a ) = -L e - 2 ^ 

Z<7i 



whereas the one for "j^ - 1 (i, 0; x\, X2) is 



u (0) (i,0;xi,x 2 ) = -^—e 
la 2 



2a 2 x 2 



As it is pointed out, the graphics included in Figure 9 correspond to the 
electric potential of the inhomogeneous case, titling the Figure with the abbre- 
viated notation (1,0) , and of the homogeneous case denoted by u\° ] (1,0). 
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Figure 11: u' 1 ' (1, 0; art, x%) 
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Figure 12: (i, 0; x\, x%) 



The graphics are traced by simply considering the value of the potentials at the 
boundary. This is, x\ — cos 8 and X2 = sin#; 6 e [0, 2w]. The same is done for 
the next Figure, where the graphics for (i, 0) and (i, 0) are presented. 

Again, as it was already expected, from a certain point of view it could be 
possible to assert that there exist a pattern on the dynamical behavior of these 
potentials at the boundary. For these cases, as the following graphics could 
suggest, the boundary potentials corresponding to the lower formal degrees, 
seem to hold the closest dynamical behavior. 



For the current vectors j W (1, 0; x%, x%) and j W (i,0;xx,X2) of (35), the 
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electric potentials are 



U 



u 



« (i, 05x1,12) 



(l,0;x 1 ,x 2 ) 




2(72X2 



1<7\X\ 



— e 




2a 1 X 1 — 1(72X2 



1 



2(72X2 



) 



The behaviors of the potentials illustrated in Figure 11 and Figure 12, could 
still suggest the existence of some dynamical relation, pointing out that it was 
already necessary to scale the graphics corresponding to the homogeneous case 
by a factor of 20, in order to show them together with such corresponding to 
the inhomogeneous case. 

Other hand, the dynamics performed in this examples can be located nearer 
of the behaviors expected when analyzing the electric potentials from a classical 
point of view: The variations of the potentials at the boundary does not hold 
any simple relation with the variations of the conductivity inside the domain 
VL. Still, by these basic examples, it is possible to suggest that perhaps from 
the point of view of the Pseudoanalytic Function Theory, applied to the study 
of the Generalized Ohm's Law, the instability of the electric potential at the 
boundary, when changes of the conductivity function inside the domain are 
taking place, could be diminished if the changes are considered for every formal 
electric potential individually. 

This could well represent an important contribution for better understanding 
inverse problems, as it is the one posed by A. P. Calderon in 1980 [3], known in 
medical imaging as Electrical Impedance Tomography. 

4 A technique for approaching separable-variables 
conductivity functions 

A very important problem is located around how to approach a separable- 
variable conductivity function once it is given a finite set of points inside a 
domain in the plane, where the electric conductivity is known (a very natural 
starting point for many physical examinations). 

This is a critical matter, since all procedures posed in this work are based 
upon the idea of possessing an explicit separable- variables conductivity function. 
Indeed, it is very difficult to find among the literature any reference about 
specific methods for interpolating separable-variables functions, even for the 
two-dimensional case. Perhaps the following idea could serve as a temporally 
departure point. 

Let us suppose we have a finite set of k% + 1 parallel lines to the xi-axe, 
dividing the domain into fci segments. Suppose we also have a set of k 2 + l lines 
parallel to the a^-axe, dividing again the segments obtained from the first step. 
The result will be a sort of grill, with a finite number k < k\ x k\ of intersections 
within the domain Q, hence we can introduce a set of points = {z' kl z'£) located 
precisely at every intersection inside J7, where z' k represents the a; 1 -coordinate 
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of Zk and z' k ' is the ^-coordinate. Finally, let us assign a conductivity value 
a (zk) to each point Zk- 

Starting with the segment of line intersecting the 22-axe at the maximum 
value reached inside f2, we can interpolate all values a (zk) corresponding to the 
points z k contained within that line. Hence we can pose a continuos function of 
the form 

fi (xi) 



at (zi) 



z" + K- 



where z'{ is the common x 2 -coordmate of all points Zk belonging to the line 
segment on which the interpolating function /1 is defined, and K is an arbitrary 
real constant such that x 2 + K ^ 0, x 2 € ft. 

The same can be done for the rest of line segments parallel to x\ found into 

O. 

We can then define the following piecewise function a (x\, x 2 ) : 



a (x 1 ,x 2 ) 



(x 2 + K) cti (xi) : x 2 E 
(x 2 + K) a 2 (xi) : x 2 E 
(x 2 + K) a 3 (xi) : x 2 E 

(x 2 + K) a kl (xi) : x 2 E 



m&xx 2 : x 2 E £1, z" — 

"+ z 2 Jl z 2+ z 



Zi - 



2 

-2 ~\~ z i 



7" — 

7" — 

z 3 




, min x 2 : x 2 E 



where z k : k 



1,2, ...,ki are the common x 2 -coordinates of every set of z k 
points belonging to the same line segment parallel to x\ inside 0. 

It is easy to see that the piecewise function a(x\ 1 x 2 ) defined in the last 
expression is a separable- variable function, so it can be employed when numeri- 
cal calculations are performed for approaching higher formal powers than those 
considered in this work. Notice also that the same idea can be extended for 
segmentations of a bounded domain fi by means of polar traces. 

Moreover, this basic idea could be also useful when considering the resulting 
electrical impedance equation when studying the monochromatic time-dependet 
case: 

div (7gradu) = 0. 

Here 7 = icue+a represents the electrical impedance, u is the wave frequency and 
e is a scalar function denoting the electrical permittivity. The reader can verify 
that virtually all mathematical methods explained in the above paragraphs can 
be extended for this complex case. 
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